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We construct models of analytic QCD (i.e., with the running coupling parameter free of Landau 
singularities) which address several problems encountered in previous analytic QCD models, among 
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I. INTRODUCTION 



In the perturbative QCD (pQCD), the running coupling apt{Q^) = cts(Q^)pt/''i': as a function of = —q^ {q being a 
typical four-momentum transfer of the considered physical process) , has so called Landau singularities at < Q'^ < 
(A^ ~ 10~^ GeV^), thus not reflecting the analyticity of the space-like observables I?(Q^) for all G C\(— cx),0] 
dictated by the causality of quantum field theories. As a consequence, the evaluated expressions of such observables 
in pQCD have wrong analyticity properties and become unreliable at low \Q'^\. In order to overcome this fundamental 
problem of pQCD, several attempts have been made to restore the correct analytic (i.e., holomorphic) properties 
of both the coupling parameter a{Q^) and the related evaluated expressions of observables, which all go under the 
generic name of analytic QCD (anQCD). 

Various models of anQCD found in the literature (some of them: Refs. for reviews and further references 

see [l]-[in]), among them the most popular minimal analytic (MA) model of Shirkov and Solovtsov [T], have faced 
criticism based mainly on one or both of the following points, one being theoretical and the other phenomenological: 

1. The analytic running coupling parameter a(Cp) differs at large from the ordinary pQCD coupling apt(Q^) = 
asiQ"^)/-^ by power terms ~ {A?/Q^Y' or - {A? / Q^Y' IvT {K^ / Q^) [in MA: by terms - [K^/Q^Y]- However, 
these then lead in (inclusive) physical observables to the corresponding power corrections which, nota bene, 
come from the ultraviolet (UV) regime [TT]. If such observables are calculated within the operator product 
expansion (OPE) mechanism, it is readily seen that such power terms are in conceptual contradiction with 
the general OPE philosophy which has been vigorously advocated in particular by the ITEP group [H] (see 
also, e.g., Ref [13]) and whose validity is strongly indicated by the success of the related QCD sum rules. This 
philosophy rests on the assumption that OPE is true in general (not only in the perturbative approach) and 
that it allows to separate short-range contributions from the long-range ones. And it is only the long-range 
contributions which should lead to power corrections, reflecting the nonperturbative physics. Thus, there is no 
space for UV-generated power corrections within the ITEP-OPE philosophy. 

2. In the widely used MA model, the prediction for one of the best-measured low-energy QCD observables, namely 
the strangeless r,-, the decay branching ratio of the r lepton into nonstrange hadrons, lies in the region r^. sa 0.14 
|14) . significantly below the experimental value ~ 0.202 ± 0.004. It appears that anQCD models in general 
tend to give too low values of r,- [T5] . 
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Point (2) (the Tr problem) can be addressed by introducing additional parameters (Refs. [ZIIH]); in MA, these 
parameters (the quark masses) have to be chosen unusually large tl6j. On the other hand, point (1) (ITEP-OPE) has 
not been addressed in a more systematic way in the literature hitherto.^ Within this letter we try to develop a version 
of analytic QCD which addresses both problems mentioned above. We base our approach on the assumption that the 
singularity structure oi a{Q^) reflects the singularity structure of space-like observables 'D{Q'^) in a "minimal" way, i.e., 
dictated by physical principles (causality and unitarity) which means that a(Q^) is analytic in the complex Q^-plane 
with the exception of a cut along the negative semiaxis starting at = —M'^^^^ < (there are no massless hadrons). 
Specifically, we expect a{Q^) to be analytic at = 0. The method of identifying such an analytic coupling consists 
in starting with an appropriate ansatz for the related beta function and reconstructing from it the coupling a{Q^) by 
solving the renormalization group equation (RGE).^ Such an approach (cf. Ref. [19]) is natural because ITEP-OPE 
condition can be implemented in this approach in a very simple way, by requiring that beta function /3(a) as a function 
of the coupling a be analytic in point a — 0, and have there the pQCD Taylor expansion /3(a) = —/3oa'^{l + cia+0{a'^)) 
where /3o and ci — /3i//3o are universal. However, having the ITEP-OPE condition easily implemented in this way, 
it turns out to be very difficult to find a beta function which gives simultaneously an analytic coupling a(Q^) (i.e., 
analytic in the complex Q^-plane with the exception of the negative semiaxis) and which gives high enough value 
Tr ~ 0.20 (i.e., compatible with the experimental measurements). In order to obtain analyticity of a{Q^), we are led 
to restrict ourselves to certain classes of beta functions. However, the obtained values of r,- turn out to be significantly 
too low unless the beta functions are further modified in a peculiar, perhaps intriguing, manner. 

In this letter, in Sec. [lljwe motivate the first class of beta functions which lead to the analyticity of a{Q^) while 
respecting the ITEP-OPE condition. In Sec. HI we modify these beta functions in such a way as to obtain the correct 
value of Tt while maintaining the analyticity and the ITEP-OPE condition. Section IV summarizes our results and 
outlines the prospects of further phenomenological applications of the obtained models. 



II. BETA FUNCTION ANSATZE FOR ITEP-OPE AND ANALYTICITY 



The renormalization group equation (RGE) 

Q^^^ = /3(«(Q^)), (1) 

determines the running coupling a{Q^) at (in general complex) once an initial condition a(Q?^) = ajn is imposed. 
We will impose the initial condition in the present anQCD versions at the scale Qf^ of the 3 — 4 active quark flavor 
threshold; we choose Qf^ = (Swc)^ (~ 14.5 GeV^). The value of ajn = a(Qfn) obtained as usual in perturbative 
QCD (pQCD). The analytic coupling we get is valid for n/ = 3 (the three active quarks u, d and s being almost 
massless), and for higher energies the standard pQCD couplings can be used because our versions of anQCD, at such 
high energies, practically merge with the pQCD due to the ITEP-OPE condition: 

Kg')-apt(Q')|<(AVQ'r (2) 

at Q2 ^ A2(- 0.1 GeV^) for aU positive n. 

More specifically, in the renormalization scheme (RSch) dictated by the expansion of our beta function /3(a) in 
powers of a (i.e., the parameters c„ = /3„//3o, for n > 2), we will require that our a{Q'^) at Q = 3mc achieves such a 
value a((3mc)^) = which leads to the value a(M|; MS) = O.IIQ/tt once we (exactly) change the RSch at Q = Sruc 
to MS and run the coupling to Q = Mz with MS perturbative RGE. The latter running is performed at four-loop level, 
taking the RGE thresholds n/ = 3i— )-n/=4at(5 = Srric and n/ = 4i— )-n/ = 5at(5 = 3mb, using the procedure of 
Ref. [5D] with three-loop threshold matching conditions (for two-loop matching conditions, cf. Refs. pTH25] ). We note 
that the value a(Mf; MS) « 0.119/7r is obtained by application of pQCD evaluations to QCD observables of higher 
energies (IQ^j > lo GeV^). 

With such a fixing of the initial condition, integration of RGE ^ in the complex Q^-planc can be made more 
transparent by introducing the new complex variable z = l'!i{Q'^ / fif^) , with /ijjj being a fixed scale; we chose = 3toc. 



^ In Ref. [3| the problem is addressed in an approximate way, by requiring the aforementioned power index n to be large (n = 4). In 
Ref. jllj it was explored whether an analytic coupling respecting ITEP-OPE can be constructed directly; it turned out to be difficult, 
and several parameters had to be introduced. 

^ In another context, an all orders beta function for non-supersymmetric Yang-Mills theories was proposed in Ref. |17| . inspired by the 
Novikov-Shifman-Vainshtein-Zakharov beta function of A/" = 1 supersymmetric gauge theories I18| . 
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The entire complex Q^-plane (the first sheet) then corresponds to the complex z stripe: — tt < Im(z) < +7r. The 
complex Q^-plane C\(— oo, 0] where a{Q^) has to be analytic corresponds to the complex z stripe — tt < Im(z) < +7r, 
while the Minkowskian semiaxis < corresponds to Imz = — vr; the point = corresponds to z = — oo, and 
— (3mc)^ to z = 0. Using the notation a{Q^) = F{z), RGE Q can be rewritten in the form 

^ = P{F{z)) , (3) 
az 

in the semi-open stripe —tt < Im(z) < +Tr, and requiring for the analyticity of a{Q^) in the sector C\(— cx),0] 
equivalently the analyticity of F{z) in the open z-stripe — vr < Im(z) < +Tr dF/dz = 0). If we write z = x + iy, 
and F = u + iv, RGE ([3| can be rewritten in term of real functions u, v and real variables x, y 

= Re/3(u + ^v) , = Im/3(. + ^v) , (4) 

= -Im/3(u + ^v) , = Re/3(. + ^v) . (5) 

dy dy 

If in ^ is an analytic function of F at F = 0, then a(Q2) fulfills ITEP-OPE condition ^. This will 

be demonstrated in the following lines by assuming that the ITEP-OPE condition is not fulfiled and showing that 
consequently /3(F) must be nonanalytic at F = 0. If ITEP-OPE condition ^ is not fulfilled, then there exists a 
positive no such that 

6aiQ^) = a(Q2) _ a^^^Q^) « k(AVQ')"° (6) 
when ^ A^. Due to asymptotic freedom at such large Q^, apt(Q^) is 

"P^^^'^ " /3oln(QVA^) +^(l^l"(QVA^)/ln^(QVA^)) , (7) 
and the power term can be written as 

(AVQ')""«exp(-if/apt(Q2)) , (8) 

where'^ K — no/Po. When we apply Q'^d/dQ'^ to relation (|6|, and use expression we obtain 

/?(a(g2))-/3pt(apt(Q2))«_„^«exp(-A7apt(Q2)) . (9) 

Now we can replace a(Q^) in the first beta function in Eq. ^ by apt{Q^) + Kexp(— if/apt(Q^)), due to relations ^ 
and ([s]), and Taylor-expand the /3-function around apt(Q^) (^ 0). This then gives 

/3(apt) -I- Kexp(-Jtr/apt) ^^^"'* | + O (exp(-2if/apt)) = /3pt(apt) - noKexp(-if/apt) . (10) 



Since this relation is valid for small values of |apt|, the derivati ve d (i(a)/da at a = apt on LHS of Eq. (10) is very small 
(about — 2/3oapt) and can be neglected. This means that Eq. (10) can be rewritten for small values oFapt — F as** 

/3(F) « /3pt(F) - mncM-K/F) . (11) 

While /3pt(F) is analytic at F = 0, the term exp(— if/F) is nonanalytic at F = 0. Therefore, the non-fulfillement of 
ITEP-OPE condition ([2| implies nonanalyticity of /3(F) at F = 0, and this concludes the demonstration. 

In addition, since at small F the beta function has to respect pQCD, the following condition must be imposed on 

it: 

P{F) = -PoF^[\ + ciF + C2F^ + ---] , (12) 



3 If the terms O (In ln(Q2/A2)/ 1^2 (^2/^2)-) ^q. ^ axe included, expression exp(— X/opt) gets replaced by 

exp(— i<'/apt)(/3o'*pt)~"''^^'''''' (l -I- 0(a\rP' a)); this does not change the argument in the text . 

Due to a typo, in the published version (JPG37,075001) the second term on the RHS of Eq. ([TT| has "-(-" sign instead of "-" sign. 
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where the parameters f3o and ci = /3i//3o are universal; at n/ = 3 we have /3o = 9/4 and ci = 16/9. 

A high precision implementation of the numerical integration of RGE ([4|-([5]), e.g., with Mathematica [21], for 



various ansatze of (3{F) function and respecting pQCD condition (12) and the ITEP-OPE condition [analyticity of 
(3{F) in = 0] then indicates that it is in general difficult to obtain a result F{z) analytic in the entire open stripe 
—77 < Im(z) < +7r. In our approach we assume that the analytic coupling a{Q^) reflects all the major analyticity 
aspects of the space-like observables ViQ^) (such as Adler function, Bjorken sum rules, etc.). This means that our 
a{Q^) is analytic even at the origin = z — — oo). This condition, in general, implies 

a(Q2) = ao + ai(QVA2) + ©[(QVA^)^] , (13) 



where Cq = a((3 = 0) = F(z = — oo) < cx^. By applying to Eq. (131 the RGE derivative Q {d/dQ ), we can see that 



the beta function /3(a) = I3{F) then has a Taylor expansion around the point oq with the first Taylor coefficient equal 
to unity 

p{F)^l^{F-ao)+0[{F~aof], (14) 

which can be equivalently expressed as 

P'{F)\f=,,=+1 . (15) 

If assuming the analyticity of a(Q^) at = in a more exceptional way a(Q^) = oq + ©[(Q^/A^)"] with n > 2, this 
implies the condition 13' {F)\p^ao — it turns out that in such cases the RGE-solution F{z) has Landau singularities, 
at Imz — -in: /n] ther efor e, we discard such a case. 



PQCD condition (12) for the universal parameters /3o and ci, the Q — analyticity condition (15), and the 



ITEP-OPE condition can then be summarized in the following form of the beta functions: 

/3(F) = -/3oF2(1 - Y)f{Y)\Y^F/a, , (16) 

where function flY) is analytic at F = (ITEP-OPE) and at y = 1 and fulfills the conditions 

f {Y) = l + + c^ao)Y + 0{Y^) , (17) 
ao/3o/(l) = 1 . (18) 



Eg. ( [l7| ) is the pQCD condition (reproduction of the universal ci), and Eq. (18) is the = analyticity condition 
( jlSj ). Under such conditions, and the aforementioned initial condition at = (3mc)^, it turns out that certain classes 
of functions /, upon RGE integration ([3|, do lead to analytic coupling F{z). Even more so, the = analyticity 
condition leads in general to solutions F\z) ~ a{Q^) which have analyticity even on a certain segment of the negative 
Q^-axis Im(z) — —n]: —M^^^. < < [o — oo < Re(z) < a;thr], Mthr being a "threshold" mass, i.e., the cut 
semiaxis in the complex Q^-plane is {—oo,—M^Yir\- 

For example, when f{Y) is a polynomial or a rational (i.e., Pade, meromorphic) function, then there exist certain 
regions of parameters of these f{Y) functions for which F{z) is analytic (o o,{Q^) analytic in the entire complex 
Q^-plane with the exception of the cut semiaxis (— oo, — Af^^j^^.]). This can be also checked and seen by analytical 
integration of RGE (|3| in such cases 

z = G{F) , G{F{z)) = / . (19) 

Jain P[-t' ) 



Namely, when f{Y) is a polynomial or rational function, integral in Eq. ( 19 ) can be performed explicitly (analytically). 
From such a solution one can see that a pole {F — oo) is attained on the negative semiaxis (at = ~-^^thr < 
i.e., aX z — xthx — ii^), and that other poles and singularitie s w ould not appear at least for certain range of values of 
the free parameters |25j. The = analyticity condition (18) turns out to be crucial for such a behavior. 



However, in this approach we encounter a serious problem: virtually all the choices of the f{Y) functions which 



fulfill the aforementioned conditions ( 17 )-( 18 1 and whose numerical solution is, at the same time, an analytic function, 
lead to too low values of the semihadronic r decay ratio (with AS" = 0): r,- < 0.16, while we need 0.20. The 

(LB) 

"leading-/3o" (LB) contribution is ri- < 0.15 for various classes of beta functions that we tried; if it is possible to 

adjust free parameters in the beta function ansatze in order to increase r)- beyond values 0.15, Landau singularities 
of the obtained F{z) [— a{Q'^)] appear. At first, for all the chosen classes of beta functions, the corrections beyond 
LB (bLB) to Tr were very small (< 0.10), and the value r^- « 0.20 could not be achieved (some elements of the 
calculation are outlined in the Appendix). 

This problem is partly a refiection of the fact that, when the analytization of the coupling eliminates the offending 
nonphysical cut < < A^ of apt(Q^), the quantity r^. tends to decrease because the aforementioned cut gave a 
positive contribution to |I5J. 
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TABLE I: Inp ut parameter values of the three considered /3-ansatze, and some resulting values of other parameters: C2, 
expansion ( 12 1, and a{Q^) at Q = irric and Q = (n/ = 3 used). 



C3 of 



/old 


input /old 


input /fact 




C2 


C3 




a ([irricf) 


a(0) 


P30 


it)i = 1 + iO.45 


;f = 43.2, B 


= 5000 


-243.6 


-250.1 


-12.00 


0.0545 


0.4596 


Pll 


Fpoic = -10. 


K = 7.0,5 = 


= 4000 


-213.3 


-293.9 


-6.44 


0.0577 


0.1995 


EE 


yi = 0.1, fci = 10.,fc2 = 11. 


if = 5.27,5 


= 1000 


-104.5 


-322.7 


-5.88 


0.0613 


0.2370 



III. R(TAU)-PROBLEM: MODIFICATION OF BETA FUNCTION ANSATZE 

Since LB contribution r\ cannot be increased further, it appears that the only way to increase the total calculated 



rr is to increase the beyond-the-leading-/3o (bLB) contributions: NLB, N LB, etc. A choice of the beta function (16) 



in our approach fixes also the coefficients C2, C3, etc., that appear in the power expansion (12) of I3{F) in powers of F . 
On the other hand, the coefficient T2 in the third term (N^LB) of the expansion of r^- beyond the LB [see Eqs. ( A14) 
and ( A17[)] c ontains a term — C2; if C2 can be made significantly negative (c2 <C —1) by a suitable choice of beta 
function ( |16[ ), then T2 and, consequently, N^LB term in expansion of will become significantly positive, increasing 
thus the evaluated value of (note that coefficient Ti of the second, NLB, term is accidentally small, Ti = 1/12, 
and independent of beta function). On the other hand, we do not want to reduce as significantly the LB contribution 
when we increase T2] and the universal ci coefficient must remain unchanged during such a modification. 
A modification which achieves the aforementioned effects is the following: 

/old(r) ^ UcAY) = /old(i")/fact(r) , (20) 

The modification factor /fact(^) is chosen in such a way (K ^ B) that, for most of the values of Y, it is close to one. 



Therefore, it does not change significantly beta function (16 1. This means that, if before the modification the LB part 



of Vt was reasonably large (say, 0.14-0.15), it will not be changed (reduced) very significantly now. PQCD condition 



(17) will not be modified by such /fact(^) because it modifies the expansion coefficients of I3{F) only at order ^ F** 
(i.e., C2) and higher. However, since 1 <^ the modification factor /fact(^) can decrease the value of C2 significantly 
and thus increase significantly the third term in the expansion of r,-. Numerical investigations indicate that this is 
really so, and that, moreover. Landau singularities are not introduced by such a modification. The latter point can 
be understood even by analytical (i.e., explicit) integration of the RGE in such a case when /om is a polynomial or a 
rational function 

The solution, however, comes at a price. The aforementioned modification increases very significantly the absolute 
values of the higher expansion coefficients c„ (n > 4) of beta function. As a consequence, coefhcients |r„| [« Cn/{n—l)] 
in the expansion become very large when n > 4. This means that the expansion series for rr starts showing signs 
of divergence after the first four terms. On the other hand, the behavior of the first four terms (including N'^LB) 
indicates reasonable behavior (similar is the behavior of asymptotically divergent perturbation series in pQCD). 

The fact that the values of parameters |c2|, Icsj, etc., are large does not mean that we are working in a "wrong" 
renormalization scheme (RSch). The specification of the RSch in terms of coefRcients Cj = /3j//3o (j ^ 2) is apparently 
a perturbative concept, applicable in the regime \Q'^\ ^ A^. It appears that our beta function /3{F) not just fixes 
a certain set of values cj {j > 2), but it reflects also certain nonperturbative aspects via its set of zeros and poles 
in the complex i^-plane. For example, the finite value ao = a{Q^ = 0) is a zero of the beta function; the function 
fncw{y) [with Y — a{Q'^)/ao — F{z)/ao] has possibly some zeros and/or poles on the real axis [but not in the 
interval Y e (0, 1)], and it has two zeros and poles on the imaginary axis close to the origin [at F — ±iB~^^^ and 
F = ±i{B + K)^^^^, respectively]. It appears that, while we might be able to go from one set of values of c/s to 
another in this framework, we cannot go to the "tame" pQCD schemes such as MS or 't Hooft RSch. For example, the 
't Hooft RSch (c2 = C3 = . . . = 0), under the assumption of the ITEP-OPE condition, gives us /3(F) = -l3oF^{l + ciF) 
and the solution in such a case violates analyticity, it has namely a Landau cut [H]). Thus it cannot be physically 
equivalent at ^ A^ to RSch's of our beta functions. The same is true for MS RSch, at least in its hitherto known 
truncated form. These considerations lead us to intriguing questions which may be clarified in the future. 



If we choose for / = /new = /old/fact in /^-function (16) for the part /om of Eq. (20) simply a polynomial, the 



LB-part of Vt- remains low unless the polynomial degree is at least three (model "P30") 



P30 : foid{Y) = (1 - wiY){l - W2Y){1 - W3Y) 



(22) 
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TABLE II: The first four terms in expansion ( A14I of r-r and tfieir sum, in tlie tiiree considered models. In parentiieses are tiie 
corresponding results for expansion (A19I. RScl parameter is C = 0; n/ = 3. The last column are variations {5) of the sums 



when RScl-parameter C increases from to ln(2). 



/o 



Id 



LB (LO) 



NLB (NLO) 



N^LB (N^LO) 



N-'^LB (N^LO) 



sum (sum) 



P30 
Pll 
EE 



0.1002 (0.0818) 
0.1065 (0.0881) 
0.1251 (0.0990) 



0.0005 (0.0100) 
0.0006 (0.0111) 
0.0007 (0.0147) 



0.0952 (0.1016) 
0.0892 (0.0961) 
0.0666 (0.0774) 



0.0060 (0.0066) 
0.0057 (0.0062) 
0.0096 (0.0107) 



0.2018 (0.2000) 
0.2020 (0.2015) 
0.2020 (0.2017) 



2.5%(2.7%) 
1.5%(1.7%) 
2.4%(2.7%) 



TABLE III: Bjorke n pola rized sum rule (BjPSR) results dsjiQ^) in the three considered models, for the sum of the first four 
terms in expansion ( A22 1. RScl parameter is C = 0; n/ = 3. In parentheses, the corresponding results for the first four terms of 
expansion ( A23 1 are given. In brackets, the corresponding variations of the results under the RScl variation are given (see the 
text for details)^ The experimentally measured values are (Ref.gT]): 0.17±0.07 for = 1 GeV^ 0.16±0.11 for = 2 GeV^; 
0.12 ± 0.05 for Q = 2.57 GeV. 



/old 


dBj(O^) : Q = 1 GeV 


Q = V2 GeV 


Q = 2.57 GeV 


P30 


0.248 (0.247) [4.8%(5.7%)] 


0.201 (0.200) [4.5%(5.3%)] 


0.145 (0.143) [3.6%(4.3%)] 


Pll 


0.218 (0.224) [2.1%(2.5%)] 


0.191 (0.194) [2.1%(2.4%)] 


0.146 (0.146) [2.8%(3.3%)] 


EE 


0.215 (0.227) [3.1%(4.2%)] 


0.188 (0.194) [2.3%(2.9%)] 


0.141 (0.141) [2.8%(3.8%)] 



For /old being a cubic polynomial, the number of free real parameters is four (two in the polynomial, and B and K in 
/fact)- This is so because initially we have six real parameters [wi, W2, w^, B, K, oq = a(0)]; two of them, e.g. and 
ao, are eliminated by the Ci condition (17) and the = analyticity condition (18). The two free parameters in the 

fLB) 

polynomial (e.g., two of the three roots) can be adjusted in such a way as to get the highest possible values of r)- ' 
(« 0.13) with /old alone (i.e., when /fact '-^ 1) while still keeping the holomorphy of F{z). Then the parameters B 
and K of /tact can be adjusted so that r-r « 0.202, the experimentally measured value. ^ These adjustments still leave 
us certain small freedom in fixing the four parameters [respecting also the condition 
the behavior of F{z) changes only little when we vary the four parameters under suet 



2T|: 1 < X < B]. However, 
conditions. In Table |lj first 

line (model P30), we present some of the results of this model for a representative choice of input parameters in this 
case: wi = 1 + iO.45 (and W2 = 1 — zO.45; as a consequence, = —3.817; wj'a being the tree inverse roots of /old); 
K — 43.2, B — 5000. In Tablcjlll first line, we present results for the first four terms of r^. expansion in the approach 
described in Appendix [Eqs. ( A12 ) and ( A14[ )] and their sum; in parentheses, the values of the corresponding first four 
terms are given in the case that no large-/3o (LB) resummation is performed. We can see that the series of r^- shows 
marginal convergence behavior when the LB-terms are resummed and three additional correction terms are included 
[see Eq. ( A14)]. If LB terms are not resummed, the convergence behavior is worse. Furthemore, the estimated value 
2.0, i.e., the series becomes divergent starting with the fifth term, 
first line, we present the results of the calculation of the BjPSR (iBj(Q^) in this P30 model for 
various values of the momentum transfer parameter Q^, taking into account the first four terms and performing 



of the fifth term is 
In Table 



III 



in parentheses, the corresponding summation of the first four terms without LB 
Eq. (A23)]. The predicted results are within the large experimental uncertainties for 



LB resummation [see Eq. (|A22 ) 
resummation is given [see 

dBjiQ^), except in the case — I GeV^ where the model predicts by about one a higher value. 

If we choose /om to be a meromorphic rational (i.e., Pade) function, it turns out that already the simplest diagonal 
Pade P[l/1] (i.e., ratio of two linear functions of Y) can do the job (model "Pll") 



Pll 



/old — 



(1 - Y/Yp,i,) 



(23) 



In this case, we have at first five real parameters {Yq, I'poio, B, K and oq], but two of them, e.g. Yq and oq are 
eliminated via the ci -condition and the = analyticity condition, Eqs. (17) and (18). We can proceed in the same 
way as in the previous case in order to (more or less) fix the three free real parameters ypoic, B, K. The results of a 



^ The value of Tr with A5 = and without mass contributions is = 0.202 ± 0.004; for details we refer to Ref. [25| : it is extracted 
from the ALEPH-measured |^-30 (V-l-A)-decay ratio iJr(AS' = 0) as in App. E of Ref. [7], by eliminating non-QCD contributions and 
the (small) quark mass effects. The result here differs slightly from the one of App. E of Ref. [7] (0.204 it 0.005) because of the slightly 
updated value of iJr(A5 = 0) = 3.479 ± 0.011 (Ref. [30]) and an updated value of |y„d| = 0.97418 ± 0.00027 (Ref. [31]). 



7 



representative choice of these input parameters are presented in the second hne (model Pll) of Tables I pTT The zero 
of fo\d{y) turns out to be at Yq = 0.6874. We see that the series for shows reasonably good convergence behavior 
in the first four terms. Inclusion of the fifth term (« —3.7) destroys the convergence, as in P30 case. Furthermore, 
BjPSR predictions now all lie within the one a uncertainties of experimental values. 

It turns out that we can choose the function /om in certain more complicated ways and fulfill all the imposed 
conditions. For example, we can choose it to be a product of exponential function of type (cxp(— y) — l)/Y and its 
inverse, both of them rescaled and translated by specific parameters (model "EE") 



EE 



(exp[-fci(r-ri)]-i) 



MY-Y2)] 



ieM-k2{Y-Y2)]-r 



X/C(fci,ri,fc2,l2) 



(24) 



where the constant /C gives just the required normalization /oid(^ = 0) = 1. At first we have seven real parameters 
(Yi, ki, Y2, k2, B, K, and Oq); two of them, e.g., Y2 and qq, are ehminated by conditions (171 and (18). We need 
< fci < /c2 to get physically acceptable behavior. It turns out that with / function being that of Eq. ( p4| (when 

f LB) 

/fact = 1, i-e., K = B = 0), the value of r)- can be increased maximally to about 0.15 while keeping F{z) analytic, 

f LB) 

if parameter Yi achieves the value Yi « 0.1. Increasing Yi further tends to increase the value of , but the 
analyticity of F{z) is destroyed through appearance of (Landau) singularities within the stripe — tt < Imz < +tt. The 

(J R) 

values of ki and ^2 have to be comparatively large and close to each other if r)- is to be kept large. Parameters B 
and K of /fact can then be adjusted so that r^. 0.202 is reproduced. 

In the third line (model EE) of Tables I pTT we present the results in this case for a representative choice of input 
parameters Yi, fci, fc2 and B and K. We see that now the convergence behavior of the series of the first four terms 
of Vj- is quite good, even when the LB terms are not resummed. Inclusion of the fifth term (« — !•) destroys the 
convergence, as in the previous two models. Furthermore, the results of BjPSR agree well with the measured results. 

In b oth T ables |TT] and |TlT we use the renormalization scale (RScl) parameter C — [cf. Eqs.( |Al4| )-( |A15[ ), ( |A19[ ), 
(A22)-(A23)]. If we vary C towards smaller values [C — ln(l/2)], the results change insignificantly, except in the case 
of BjPSR at g2 = 1 GeV^ in Pll and EE. If we increase C to ln(2), the results decrease, and the percentages of such 
decrease of are given in the last column ("(5") of Table |ll) and for BjPSR dBj(Q^) are given in Table III in brackets. 
Only in the case of BjPSR at = 1 GeV^ in Pll and EE these percentages mean the variation (decrease) of the 
result when C goes down to ln(l/2). In parentheses, the corresponding values are given when the LB terms are not 



resummed, cf. Eqs.(A19l and (A23|. If only three terms are included in our calculations, the variations of the results 



for Tr and BjPSR under the aforementioned variations of RScl significantly increase, in general to about 10%. 

If we use as the basis of our calculations of and c?Bj the truncated expansions in powers a" like Eq. ( A3 ) , instead 
of the truncated expansions ( A4 1 in logarithmic derivatives a„ ( A5 ) , the results turn out to be significantly more 
unstable under the variation of RScl. E.g., the value S in Table |Il| in the case Pll changes from 1.5%(1.7%) to 
8.2%(9.9%), and the value of changes from 0.2020 ± 0.0031 (0.2015 ± 0.0034) to 0.2815 ± 0.0232 (0.2747 ± 0.0272). 




s 



10 



12 



14 



-a.05 



0.0 



0.1 



FIG. 1: (a) Analytic coupling a(Q'^) (full line) at positive < ' 
<12(Q^) (dashed line) and as (Q^) (dot-dashed line) [cf. Eqs. \ 
as in (a), but at very low > 0. 



0.2 0.3 



0.4 



O.S 



< (3mc)^, in model EE; included are also higher order analytic couplings 
lSI], for better visibility scaled by factors 5 and , respectively; (b) same 



All these results show that model EE is very similar to model Pll, but significantly different from model P30. 
Further, the threshold values Xthr in models EE and Pll are similar (see Table [T]): Xthr ~ —6.; this corresponds to 
the threshold mass Q^j^j. = — M^^j^^. for the discontinuity function pi(cr) with values Mthr = (3toc) exp(xthi /2) ~ 0.2 




EE: Y 1^.1, kl^lO., k2= 1 1, K=5.27, B =10^3, nfrJ 

v[f:.-!r] 
4r 



(a) 



EE: Y 1^.1, kl^lO., k2= 1 1, K=5.27, B =10^3, nfrJ 

n[\.-!f] 



(b) 



-S -6 -4-2 2 

FIG. 3: (a) Imaginary part lmF{z = x — in) = v{x, 



ln[c/(3 n.icj''2] ' — ' — ' — ^ 

-S -6 

—it) of the analytic coupling F{z) 







ln[c/(3 n.icj''2] 



v{x,—it) = Ima(Q^ = —a — it) = pi(o') is the usual discontinuity function of the analytic coupling, where a 
rric = 1.27 GeV). (b) Same for the real part ReF(2 = x — in) = u(x, —n). 



a{Q^) in model EE, as a function of x. Here, 
^exp(x) (/lin = 3mc, 



GeV. On the other hand, m model P30, Xthr is much more negative: a;thr ~ —12-, corresponding to Mthr ~ 0.01 GeV. 
For all these reasons, we will consider models EE and Pll as two viable models of analytic QCD which fulfill the 

conditions imposed at the outset of this letter. 

In Figs.flla), (b) we present a{Q^) and the higher order couplings aj{Q^) {j — 2,3) [cf. Eq. (A5)], in model EE, as 
functions of at low positive < ^f^. The Figure indicates strong hierarchy a{Q^) >• 02 (v^) ^.^{Q^) >• . . . at 
all positive values of Q^. In Fig.[2]we present the three-dimensional image of \l3{F{x + iy))\ as a function of x and y; we 
can see that there are no singularities of this function inside the z stripe — tt < y(= Im(2;)) < +n; the only singularity 
is at the threshold value zthr = —5.8754 — in which corresponds to « —(0.202)^ GeV on the negative Q^-axis. In 
Figs. [3] (a), (b) we present the behavior of the imaginary and real part of the coupling F on the edge z — x — in [i.e., 
on the negative axis: — —fj,?^exp{x)]. We see the threshold-type behavior at zthr — a^thr — in = —5.8754 — in. 
The fact that these latter curves have no (step-like) discontinuities at a; ^ Xthi- is an additional numerical indication 
that the function F{z) has no singularities within the stripe —n < Imz < n, i.e., no Landau singularities. 



IV. SUMMARY 



We investigated whether it is possible to construct analytic versions of QCD which obey the ITEP-OPE principle 
of no UV-contributions to power term corrections to pQCD (A^/Q^)" and, at the same time, do not contradict the 
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measured value of the semihadronic t decay ratio Tt (which is by far the most precisely measured low energy QCD 
quantity). We constructed such models by choosing specific forms for the RGE beta- function, and found that the 
answer is positive: such theories do exist. However, the obtained solutions came at a price, because the obtained series 
for rr show divergent behavior starting with the fifth term of the series. This was so because we had to introduce 
poles and zeros of the beta function on the imaginary axis relatively close to the origin (in the complex plane of 
the coupling), in order to increase the value of r^- One model contained a cubic polynomial, another a simple Pade 
P[l/1] function, and yet another model a combination of exponential functions of the type (exp(— y) — 1)/Y . The 
last two models show better apparent convergence behavior of r^- (in the first four terms) and agree well with the 
(less precisely) measured values of the Bjorken polarized sum rule at low energies. The last two models appear to 
be numerically very similar to each other. We intend to use these two models in the future evaluations of various 
physical quantities with the OPE approach. This approach can be applied with the presented analytic QCD models 
since the latter respect the ITEP-OPE philosophy. For example, higher-twist contributions to the Bjorken polarized 
sume rule may be substantial. Such contributions were ignored in the numerical analysis here, but should eventually 
be included. 
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Appendix A: Expansions and resummations of observables in analytic QCD 



Here we refer to and summarize the approach decribed in our previous work [7j . The massless stangeless ( A5 — 0) 
semihadronic r decay ratio can be expressed in terms of the current-current correlation function H((5^) (massless, 
V-V or A- A) as 

rr^- r 4 fi - 4 V f 1 + i"in(g2 = -s) . (Al) 

TT ./n m± \ rat / \ rat j 



This integral can be transformed, via the use of Cauchy theorem in the Q^-plane^ and the subsequent integration by 
parts, to the contour integral [3H 



1 

2^ 



dcj, (1 + e^'^Yil - e'^) dAdi(Q' = ra^'^) , (A2) 

where (iAdi(Q^) = —dIl{Q^)/dlnQ^ is the (massless) Adler function whose perturbation expansion is 

00 

dAdi.(Q') = a + ^d„a"+i (A3) 

00 

= a + ^ d„a„+i . (A4) 

n=l 

Here, the coupling parameter a — a{p'^; C2, C3, . . .) is at a chosen RScl /i^ and in a chosen RSch (c2, C3, . . .) (c„ = /3„//3o), 
as are the coeffficients dn and dn' d„ = (i„(C; C2, . . . , c„_i), dn = (i,i(C; C2, . . . , c„_i). Here, C is the dimensionless 

RScl parameter: C = ln(^^/Q^). 

The higher order couplings a„+i appearing in (A4| are 

«„.r(.^)^y,^^, (n = l,2,3,...). (A5) 



® In perturbative QCD (pQCD) this use of Cauchy to relat ion l |Al|l is formally not allowed, due to the unp hysic al (Landau) cut of 
npt(Q^) along the positive axis 

< Q2 < a2; in pQCD, |Ai} and ( |A2| l are in principle two different quantities, (jA2| being the preferred 

one. 
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The two expansions in ( A3 1 and ( A4 1 are in principle equivalent (not equivalent in practice, when truncation used) , 
because of the relations 



02 = a"^ + aa^ + C2a^ + 0{a^) 

as = + {5/2)cia* + O(a^) , 04 = a"* + 0{a^) , etc. 



and the consequent relations between c?„ and dm's 



di = di , ^2 = ^2 - cidi , 

^3 = ^3 - (5/2)cid2 + [(5/2)c2 - C2] di , etc. 



(A6) 
(A7) 



(A8) 
(A9) 



The leading-/?o contribution (LB, in Refs. [3 [7] named leading-skeleton LS) to the massless nonstrange ratio was 
given in Ref. |7] in Appendix C, Eqs. (C8)-(C11 ), us ing results of Refs. [34j[35]. It is the contour integration (A2) of 

the LB-part d^^^^ of Adler function expansion (A4). While the LB part was written in Refs. [6l [7] in terms of the 
Minkowskian coupling 2li 



.(LB) 



dt 



(AlO) 



where C = —5/3, the characteristic function F^{t) is given in Eqs. (ClO)-(Cll) there, ^ and the Minkowskian (time- 
like) coupling 2ti((T) is related with the discontinuity (cut) function Pi((t) of the coupling parameter a [pi((t) = 
Ima((5^ = —a — it)] in the following way: 



-^2li(a) = --pi((7) 
diner TT 



(All) 



Since the discontinuity function is Pi{(y) = Im_F'(z) for z = \'n-{a / nf^) — it:, it is obtained as a direct byproduct of the 
integration of RGE ([3|. Thefore, it is convenient to express LB contribution (AlO) in terms of pi{a) instead of 2li(cr). 
This can be obtained from relation (AlO I by integration by parts and using relation (All I 



.(LB) _ 



dt 



Fr{t) pi{te^ml) , 



where 



Frit) 



* dt' 



jFrit') 



(A12) 



(A13) 



bmce F^ (f) consists of powers of t' and polylogarithmic functions of t' and 1/t', it turns out that integration in 
(A13) can be performed analytically. Explicit expressio n for Frjt) will be given in Ref. [5S]. Here we only mention 
that Fr{t) 1 when t — > -l-oo, and that integration in (A12) starts at a positive tthi = (Mj^jjj./m^) exp(— C), due to 
the threshold behavior of pi{(y) in our presented models. 

A systematic expansion of beyond the LB can then be written as r-r plus contour integrals of a„+i's (n > 1) 

00 

= r (LB) + J2 T„/(5„+i , C) , (AM) 

n=l 

where 



/(a„+i,C) 



1 

2tt 



(l + e'*)'(l-e'^)S„+i(e'-my^) 



(A15) 



A typo appears in the last lino of Eq. (Cll) of Ref. in a parenthesis there instead of a term -|-3 should be written +3t-^; nonetheless, 
the correct expression was used in calculations there. 
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C is an (arbitrary) renormalization scale (RScl) parameter (|C| ^ 1), and the coefficients Tj are 



T2 



Ti 



.(1) 
-10 



1 

12 



T2 + 2/3oC c^J,^ 



(C2 - C2) 



+(C2 - C2) ((5/2)ci - - 3^o(4i^ + C)) - (l/2)(c3 - C3) 



(2)^ 



(A16) 
(A17) 

(A18) 



The overHnes indicate the corresponding quantities which appear in the MS RSch with RScl parameter C = 0; 
coefficients c^J-* are determined by the /Jg-expansions of the perturbation coefficients of the massless Adler function 



d{Q'^); for details see Ref. [7], particularly Appendix A.* In particular, for Uf = 3: 



,(1) 



.(2) 
-10 



= -278.673, c 



(2) 
11 



= 1/12, 



„(i) 



0.691772; 



59.2824. The N^LB coefficients T3 and can now be calculated exactly because the N'^LO 
perturbative coefficient of the massless Adler function is now known exactly [36]. In our case {nj — 3) it turns out 
that T 2 = - 12.2554 and T3 = 1.55291. 



Eq. (A17| indicates that N'^LB coefficient T2 becomes large positive [and thus the N^^LB term in expansion (A14) 
becomes significant positive] if the beta-coefficient C2 becomes negative: C2 <C —1. Futhermore, if |c4| is large and 
dominant (as it is in our models), Eqs. (A16|-(A18) indicate that T4 « — (l/3)c4 and thus IT4I is large 



If no LB-resummation is performed in 7v(<=> in dAdi.); then r-r is obtained by performing contour integration (A2) 
term- by-term for the sum ( A4 ) 



r^ = /(a,C) + ^ d„/(a„+i,C) 



(A19) 



In practice, we have to truncate sums (A14) and (A19), by including rtmax = 3 because only the first three coefficients 
dn dn) are known exactly 



Bjorken polarized sum rule (BjPSR) dBj(Q^) is yet another QCD observable with measured values (although much 
less precisely than r^) at low energies. It can be calculated in a similar way. Its perturbation expansion can be 
organized in two ways, like in Eqs. (|A3[) and (A4) for the Adler function. LB-resummation 



2n(LB) 



dt 



FBi{t)a{te^Q^) , 



can be performed with the characteristic function obtained in Refs. [7] 



FBiir) 



W) r < 1 



r 1 



Inclusion of terms beyond the LB gives 

00 

dBj(Q') - dBj(Q')^^^^ + ^(TBj)„S„+i(e'^Q2) , 



(A20) 



(A21) 



(A22) 



where coefficients (T'Bj)„ are analogous to coefficients r„ of Eqs. (A16)-(A18), but this time based on the BjPSR 
perturbation coefficients (dBj)fc {k = 1, . . . ,n) instead of dk of Adler function. The perturbation coefficients (rfBj)i 
and (dBj)2 are known exactly [3H|, and for ((iBj)3 we use an estimate given in Ref. [JD] for Uf = 3: (dBj)3 ~ 130. 
If LB resummation is not performed, the resulting expression is 



dBj(Q') = a{e^Q^) + ^(dBj)„a„+i(e'^g2) , 



(A23) 



In Ref. [7], notation An was used instead of a„, and tn+l instead of T„. The power analogs An constructed in Refs. [6] [7] reduce 
to powers a" here because /3(a) here is analytic in a = (as a consequence of ITEP-OPE condition). In Eq. (A18) of Ref. [7] there 
is a typo, in the first line the last term there should be —Sb2i3{c^j^i + C) instead of —5b2i3cfj^i . The correct formula was used in the 
calculations there; e.g., Eqs. (89)-(92) in Ref. [7], which follow from Eq. (A18) there, are correct. 



12 



where the perturbation coefScients (rfBj)n are evaluated at the chosen RScl /i^ = exp(C)Q^ and in the RSch (c2, C3, . . .) 
dictated by /3-functions of our analytic QCD models. 
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